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ABSTRACT 

We investigate the hydrodynamics of accretion channehed by a dipolar magnetic field 
(funnel flows). We consider situations in which the electrons and ions in the flow 
cannot maintain thermal equilibrium (two-temperature effects) due to strong radia- 
tive loss, and determine the effects on the keV X-ray properties of the systems. We 
apply this model to investigate the accretion shocks of white dwarfs in magnetic cat- 
aclysmic variables. We have found that the incorporation of two-temperature effects 
could harden the keV X-rays. Also, the dipolar model yields harder X-ray spectra 
than the standard planar model if white dwarf is sufficiently massive (> IMq). When 
fitting observed keV X-ray spectra of magnetic cataclysmic variables, the inclusion of 
two-temperature hydrodynamics and a dipolar accretion geometry lowers estimates 
for white-dwarf masses when compared with masses inferred from models excluding 
these effects. We find mass reductions < 9% in the most massive cases. 

Key words: accretion, accretion discs — hydrodynamics — shock waves — stars: 
binaries: close — stars: white dwarfs — X-rays: binaries 



1 INTRODUCTION 



Field-channelled accretion occurs in a vari ety of stellar 

objects (e.g. iKonigll Il99ll: 
g et al.l 200d: ^Lamzin et al.l 



systems, from youn g stellar 



Hartman et"!!? 1993; 'Gullbring 
20011: iRo manova ct al. 2003 



Stelzer fc Schmi t^ |2004 



Gregory et al.l 120061) to compact stars (e.g. lElsner fc Lamb 
19771; iGhosh fc LambI [l978l: lAronj 1 199^; iLovelace et al 



Koldoba et al 



1995l:lLi et al.lll996l:lLi fc Wilson' 1999: Kr yukov et aLll2000l: 



2OO2I : IXoropina et al.i i2003l : ICanalle et al l 



20051 ). where the magnetic stress of the accretor exceeds the 



region decelerate abruptly to attain subsonic speeds. Typ- 
ically, the shock temperature is T ~ SGM^mn/SkR^ ~ 
10 — 50 keV (where G is the gravitational constant, k is the 
Boltzmann constant, mn is the hydrogen atomic mass, Mw 
is the white-dwarf mass, _Rw is the white-dwarf radius). The 
shock-heated flow is cooled by emission of free-free and line 
X-rays, and optical/IR cyclotron radiation. 

The post-shock regions of white dwarfs in mCVs are 
stratified in temperature and density, and the flow struc- 
tures axe in_general_w^l^^ 



ram pressure of the accretion flow. For accretion on to white 
dwarfs in magnetic cataclysmic variables (mCVs) (|Warneij 
I1995I : ICropperlll990l ) , the magnetic field not only directs the 
flow but could also dictate the radiative-loss processes (see 
e.g. IWu et al.. 2003). This, in turn, alters the flow hydro- 
dynamics, leading to observable consequences, for example, 
in the optical/lR polarisation and keV X-ray spectra of the 
system. 

The keV X-ray emitting regions of mCVs are usually 
located at the magnetic fleld footpoints on the white-dwarf 
surface. Here an optically thin plasma slab heated by a 
strong shock, at which the supersonic flows in the upstream 



(lAizu|l973l:IChevalier fc Imamurall 19821 : IWulll994l : IWu et all 
11994 ICropper et al.lll999l : IWull2000l ).These models assume 
that the electron and ions are in thermal equilibrium, via 
collisions, and share the same local temperature. (Here- 
after, these models are denoted as standard one-temperature 
(IT) models.) The standard IT models have successfully ex- 
plained the line features in the keV X-ray spectra of mCVs. 
However using the standard IT models to generate spectral 
fits for mCVs tends to require white-dwarf masses systemat- 
icall y higher than the masses measure d by other techniques 
fsee lRamsay et al.lll998l : lRamsay|2000l ). This has led to ques- 
tions on the validity the IT approximation. 

It has long been recognised (jFabian et al.l Il976l : 
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iKing fc Lasotalll979l : iLamb fc Masters! 1 19791 ) that the elec- 
trons and ions in the post-shock flows ol mCVs might be 
unable to maintain thermal cquihbriu m solely via electron 
ion collisions (see also Ima mura 1 981i; Imamura et alj 



19961: IWoelk fc Beuermannll 19961 : ISa^ton fc Wull 19991 



1987 



2001 



Fischer fc BeuermannI 1200 ll ). The departure of the elec- 



tron temperature from the ion temperature can be se- 
vere if the radiative-loss timescale of the flow is much 
shorter than the dynamical timescale. In such situations, 
two-temperature (2T) effects on the X-ray spectra are no 
longer negligible. Two-temperature effects on X-ray spec- 
tra of mCVs wer e previously invest igated assuming pla- 
nar stratifled flow jSaxton et al.ll2005l ). However, the accre- 
tion flows in mCVs are not planar. There is evidence that 
the accretion flows form curve d funnel /curta i n-like columns 
along the magnetic field lines (|Ferrariclll996l : iHeerlein et al.l 
1 19991 ). The hydrodynamics in funnel flow and in planar are 
qualitatively different, especially for tall post-shock columns 
where a variation of the gravitational potential is significant 
l|Canalle et al.ll2005l ). 

Here we investigate 2T accretion in mCVs where the 
flows are guided by the white-dwarf magnetic fleld geome- 
try. This work is a generalisa t ion of both the 2T planar flow 
calculations of ISaxton et al.l (|2005|) and the IT dip ole-fleld 
channelled flow calculations of Canalle et al.l (|2005l ). We or- 



ganise the paper as follows. In §2 we present the model as- 
sumptions and the hydrodynamic formulation; in §3 we show 
the hydrodynamical structures and X-ray spectra obtained 
from the calculations and discuss the properties of the flows 
(and their dependence on system paramters). In §4 we sum- 
marise the major findings. 



2 HYDRODYNAMIC FORMULATION 

For a gas with adiabatic index 7, density p, velocity v and 
pressure P, subject to a gravitational field g, the general hy- 
drodynamic equations for conservation of mass, momentum 
and energy can be expressed as 



^pv -I- V ■ pvv + VP 
and 



Pg 



(1) 



(2) 



i- 



V- V P 



tP /9_ 

p \dt 



V- V 



-(7-l)A, (3) 



where A is the volumetric energy loss rate due to radia- 
tive cooling. Emi ssion by electrons fa r exceeds emission by 
the ions. Unlike ICanalle et al.1 (|2005l ). who assumed tem- 
perature equilibrium between electrons and ions, we will 
consider cases where electron radiative cooling is significant 
compared to coUisional heating by the ions. This entails a 
separate energy equation for the electron sub-fluid, 

(|+v.v)p.-2^(| + v.v)pM.-i)(r-A),(4) 

where Po is the electron partial pressure. The electron 
cooling function, A = A(p, Pc), and the coUisional heat- 
ing function, V = V{p,P,Pe) , depend on local hydrody- 
namic variables llSpitzeil Il962l : iRvbicki fc LightmanI 1 19791: 
llmamura et al.lll996l ). 



We consider accretion flows that are channelled by a 
dipolar magnetic field centred on the accreting stellar object. 
The coordinate w measures paths along a magnetic field line 
towards the accretor. A transverse coordinate u = siv? Ot, is 
related to the colatitude of the accretion hot-spot {9»). The 
azimuthal angle is the same as in cylindrical coordinates. 
We seek stationary solutions for the flow structure, so we 
omit temporal derivatives from I!}-©. We recast the equa- 
tions in curvilinear coordinates, {u,w,ip). For flows along 
the field line, there is no velocity component in the u or ip 
directions. Then the hydrodynamic equations for mass flux, 
momentum flux, total energy and electron energy simplify 
to 



_d_ 

dw 



hihspv = 



hihspv = C , 



V dv 1 dP 
h2 dw h2p dw 



V dP 
h2 dw 

V dPo 
h2 dw 



7 Pv dp 
/i2 p dw 

7 Pgv dp 
h2 p dw 



-(7-l)A 



= -(7-l)(A-r) 



(5) 
(6) 
(7) 
(8) 



where C is a constant proportional to the mass flux onto 
the stellar surface, C = mfti«/i3t. The functions hi{u,w) 
and /i3(w, w) are metric terms of the curvilinear coordinate 
system (see Appendix [XJ . The term = g • i& is gravita- 
tional acceleration along the magnetic field line (i.e. in the 

w direction). 

As expressed in lSaxton et al.l (|2005| ). the rate of energy 
exchange from ions to electrons has the form 



r = Xp^/'p-^''^(p-xPc) 



(9) 



where X and x s-re numerical constants that depend on the 
plasma composition. The volumetric cooling rate of gas due 
to thermal bremsstrahlung radiation is approximated as 



Ab 



Ap^ 



(10) 



with A being anoth er constant depending on composition. 
ISaxton et a"lT(|2005l ) gave derivations of A, X and x and their 
values for a plasma of solar metallicity. We assume a total 
cooling function that includes bremsstrahlung and cyclotron 
contributions. 



A = Abr + Acy = Abr [1 + es,fcy(p, Pc, 



(11) 



with a cyclotron /bremsstrah lung emiss iv ity ratio g i ven b y 



the function f ey defi ned by IW 
Saxton et"aD ( 19971) , refined by 



1994h: iWu et al 
Saxton fc Wul 



1994); 



19991 ): 



Cropper et al] (|l999l ): ISaxton et al.1 (I2005D and extended to 
the dipolar accretion model in ICanaUe et all (|2005| ). 

We rewrite the flow density, total pressure and electron 
partial pressure in terms of the h functions, the flow constant 
C and variables that have the dimensions of velocity. 



C 



hihsv 

C 
h\h-i 



Pc = 



(12) 
(13) 
(14) 
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where ^ me asures the ratio of m omentum to mass fluxes (in- 
troduced in lCropper et al ]|l999l ) and a = Pe/Pi, the electron 
to ion partial pressure ratio. For the sake of further brevity, 
we deflne 



C 



Po = 



cr + 1 



(15) 



We separate the finite and terms in the electron-ion 
energy exchange function and the bremsstrahlung radiative 
cooling functions, 



X 



c 

hihs 
C 



a + 1 — 



C 



ap 



l/2,,5/2 



3/2 



h\h2hz 
C 



-Tv 



All 



-5/2 



-3/2 



(16) 



(17) 



/!.l/l2^3 

The functions f = hrhih^Vv'^/'^ /C and A = hxh^h^Kv^l'^ jC 
lack explic it v dependencies an d are finite everywhere. 

As in ICanalle et all ()2005l ). the equations for momen- 
tum and total energy provide 

dw V 



dA^ _ - [(7 ^ l)At;-3/^ + 7H(e - v)v + h^gu,] 
dw 7(^ — v) — V 

Definition (|14p . equations (O and (|19|) yield 



dPc 
dw 

da 
dw 



^ dw 
(7-l)(a 



dv \ 
dw J 



HPc + 



(P - Pcf da 



dw 



(C - v)vy 



±i) [.A-(a+l)f] 



(18) 
(19) 

(20) 
(21) 



Strong shock conditions apply at the outer boundary, 
while the condition v — applies at the stellar surface 
(r = 1). Our particular method of numerical solution is 
described in Appen dix [B| W e adop t the white dwarf mass- 
radius relation of iNauenberd l)l972h . 



3 RESULTS 

3.1 Structure of the post-shock flow 

Table [1] summarises the basic properties of several models 
of mCVs, with nearly direct (approximately radial) accre- 
tion onto a magnetic pole (colatitude ^* — 0.001°) and hot- 
spot area 10^^ cm^. We consider white dwarfs in the mass 
range 0.7 — I.2M0, with surface magnetic field strengths 
Bt, ^ 50 MG. The ratio of electron to ion temperatures at 
th e shock, (re/Ti) a, is t r eated as a free paranieter (e.g. as 
in Imamura et al.l 1 19961 : ISaxton fc "Wulll999l : ISaxton et al.l 
I2OO5I ). This ratio has little effect on the shock height, 7?s. 



This implies that the shock position is insensitive to the de- 
tails of heating at the shock precursor, but depends mainly 
on the white-dwarf's properties, the infiow rate and field 
geometry. 

Figure [1] shows how the shock radius varies with the 
accretion colatitude 6*,, the magnetic field strength and ac- 
cretion rate. These illustrative models assume that To — Ti 
at the shock, and accretion hot-spot a rea of 10^^ cm^ . In all 
calculations (as in lCanalle et al.ll2005l ) the shock is closer to 
the stellar surface when rh is greater. The greatest variation 
of _Rs with 8, occurs in systems where the shock radius is 



Table 1. Parameters and properties of a set of representative 
accretion shocks, varying the electron to ion temperature ratio at 
the shock (To^s/T; g). The white-dwarf surface magnetic field is 
given by Br = S«/10 MG, and the specific accretion rate rh is in 
units of g cm~'^ s~^. If the magnetic field at the hot-spot is B» 
then its magnitude at the pole is Bp = 2iJ* / \/4 — 3u. 
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large compared to R^. We find that the IT model (dotted 
curves) predicts lower shock positions than the 2T model 
does. The difference is proportionally greater in large-B* 
cases, where the shock is typically nearer the stellar surface. 
As the 2T model predicts increased Rs values, the shock 
temperatures are lower, Ts cx R^^- This effect is greatest for 
large Mw. 

Figure[2]portrays the radial distribution of electron den- 
sity (rio) and temperature (fcTg) in accretion models with 
6** = 0°, Mw = I.OM0 and m = 2 g cm"^ s"\ The IT 
and 2T models predict similar density and thermal struc- 
tures in the lower field cases (B* — MG). However for 
stronger magnetic fields {B, = 50 MG) the models predict 
considerably different shock heights and qualitatively differ- 
ent density and temperature distributions. In 2T structures, 
the hottest electrons are confined to a narrower region near 
the shock than in IT shocks with equivalent system parame- 
ters. Plotting the Tc versus rio distributions, as in the lower 
panels of Figure O shows that in 2T models with strong 
magnetic fields the density of the hottest gas is lower than 
in IT models. This will have consequences for the post-shock 
X-ray emissions. 



3.2 X-ray emissions 

For each of the solutions for the post-shock accretion flow, 
we calculate an X-ra y spe ctrum followi n g the methods of 
ICropper et al] l| 19991 ) and ISaxton et akl (|2005l '). The post- 
shock volume is divided into > 10* strata according to the 
spatial steps of the numerically integrated flow profile. The 
volume of each stratum is inferred from the step size dw, 
the hot-spot area and the crossectional area scaling fac- 
tor hihs/hi,hsf The local electron density and tempera- 
ture (rie. To) are obtained from the dimensionless flow model 
and the scaling units, (i?w,Vw,C). We assume that the X- 
ray line and continuum emission are optically thin, and in- 
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Figure 1. Effect of the accretion colatitude (S*) upon tlie sfiock radius, Ts = Rs/Rm in units of the stellar radius, for cases with different 
specific accretion rates (m), white-dwarf mass (Mw) and magnetic field (S*). The accretion hot-spot area is lO^'' cm^. Solid Unes mark 
the results of the 2T model; dotted lines denote the corresponding results of the IT model. Within each panel, paired curves from bottom 
to top represent cases with field strengths B* = 50, 30, 10 MG; a fourth pair shows cases with neglible cyclotron cooling (yet channelled 
by a dipole funnel). These low-field cases are insoluble when rii = 1 g cm~^ s~^. and M„ = 1.0, 1.2Mq. 



voke the XSPE C implementation of the MEKAL therma l 
plasma model l|Mewe et al.1 1 19851 : iKaastra fc Mewd \l99^ ) 
to calculate emissivities for local {ne,Te) values. Volumet- 
ric integration over all the post-shock strata yields a syn- 
thetic X-ray spectrum. We exclude the densest strata, with 
rie 10^* cm~^, which the spectral model cannot handle 
(and which in any case represent optically thick layers merg- 
ing into the atmosphere of the white dwarf). We omit the 
effect of shadowing by a white dwarf which obscures its own 
accretion inflow. 

ISaxton et al.1 l|2005l ) studied 2T effects in a planar accre- 
tion model, without the effects of gravity and varying width 
of the magnetic accretion funnel. They found that 2T effects 



change the keV X-ray spectra most significantly for cases 
with larger M„. For given system parameters, the 2T model 
gave a harder continuum but relatively weaker line emission. 
The ratio of electron to ion temperatures at the shock was 
found to have little effect on the spectra. This was explained 
in terms of the concentration of X-ray emission near the 
stellar surface (well downstream from the shock) by which 
point coUisional energy exchanges have nearly equilibrated 
the electrons and ions. 

Our present calculations, with a dipolar accretion fun- 
nel and 2T effects, typically produce spectra such as those in 
Figure [3] The slope of the continuum is steep, but there are 
variations in the curvature of the continuum and the details 
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Figure 2. Profiles of electron temperature (upper row) and electron number density (middle row) for accretion onto 1.0 Mq white 
dwarfs, as a function of height above the surface. The left and right columns respectively show cases with negligible cyclotron cooling, 
and with dominant cyclotron cooling. The accretion spot is at the pole (6» =0°). Solid lines mark the results of the 2T model; dotted 
lines denote the results of the IT model. Electron and ion temperatures match at the shock and stellar surface, but the electrons are 
cooler elsewhere. The bottom row shows the temperature/density phase structure: in high-field cases the 2T model has lower densities 
in the hotter strata. 



of line strengths. The effects of the various system param- 
eters are best iUustrated by plots of the ratios of spectra 
calculated with different conditions or assumptions. 



3.2.1 Dipolar vs planar accretion 

Figure |4] shows the ratios of X-ray spectra calculated in 2T 
models with identical system parameters but dipolar verses 
planar accretion geometry. In all the cases we calculated, the 
dipolar model provides a harder continuum, and reduced line 
emission, particularly for photon energies < f keV. These 
effects are small for a white dwarf of O.JMq, but are more 



considerable for I.OA/q and more massive cases. The dipo- 
lar (2Td) and planar (2Tp) models are more in agreement 
regarding continuum at photon energies < 1 keV if the ac- 
cretion hot-spot is at high S« (further from the magnetic 
pole). 



3.2.2 Two-temperature vs one-temperature dipolar inflows 

Figure [5] shows the quotients of X-ray spectra over the 
0.2—10 keV band, from IT and 2T dipolar calculations for a 
white-dwarf mass of O.TM©. The X-ray continuum is gener- 
ally brighter in a 2T model. In cases where bremsstrahlung 
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is the dominant cooling process (B* = MG) the X-ray 
continuum has equivalent spectral slopes in ITd and 2Td 
models, up to photon energies ~ 5 keV. Above this energy, 
the ITd model predicts a marginally harder continuum. The 
ITd models predict ~ 0.5% stronger emission lines around 
1 keV, but lines around ~ 7 keV are weakened by a similar 
degree. 

For non-zero but small magnetic fields (_B» = 10 MG) 
the shape of the continuum is similar to the zero-field cases, 
but the ITd and 2Td calculations yield greater discrepancies 
about the line strengths, up to ~ 1%. The specific accretion 
rate, rh, has a slightly greater efi'ect on the disagreements 
between the predicted lines. 

For a stronger magnetic field [B, = 30 MG) the ITd 
model yields a softer continuum than equivalent 2Td calcu- 
lations. The ITd model generally predicts stronger emis- 
sion lines, by up to ~ 2% for the cases with lowest m. 
Two-temperature effects are less apparent in the keV X-ray 
continuum and lines when the accretion rate is larger. The 
spectral effects generally strengthen as the magnetic field 
increases; the greatest extreme in our plots is the case with 
B, = 50 MG and m = 1 g cm"^ s"\ 

The discrepancy between spectra predicted in the ITd 
and 2Td models is more acute for systems with larger white- 
dwarf mass. For cases with Afw = l.OM© (see Figure |6]) the 
oversoftness of the ITd predicted continuum is more pro- 
nounced in the B* = 30 MG case than in the comparable 
B, = 50 MG cases with Mw = 0.7Mq. The lTd/2Td pat- 
tern of over- and under-prediction of line emission between 
~ 7 keV and ~ 9 keV varies considerably with B« and rh 
(see the upper right panels of Figure |6]). 

For systems with a white dwarf with Afw = 1.2Af0, the 
discrepancies are even greater than for Afw = l.OAf0 cases: 
the continuum around ~ 10 keV can be around 10% under- 
predicted in the ITd model, relative to the 2Td calculation 
(Figure [7]). A downward spike in the H-like Fe Ka emission 
that appears for l.OAf© cases (in Figure |6l representing un- 
derprediction of this line by the ITd model) is absent when 
Afw = 1.2Af0, with all other parameters equal. 

3.2.3 Latitude dependence 

The X-ray spectrum also varies with the location of the ac- 
cretion spot on the white-dwarf surface. Latitudinal efi'ects 
are slight when the shock height is small compared to the 
stellar radius (e.g. Af„ — 0.7Mq, m = 5 g cm~^ then 
we obtain < 2% enhancements in the ~ 1 keV lines. For 
the same white-dwarf mass but lower accretion rate (e.g. 
Afw = O.7Af0, m = 1 g cm~^ a low-latitude accre- 

tion stream provides a softer X-ray continuum than polar 
accretion. Figure [S] shows the ratios of spectra: for inclined 
accretion cases (0, = 10°, 30° and 60°) compared to vertical 
accretion onto the pole (S« =0°). For low-latitude accretion, 
lines are enhanced up to < 20% at photon energies < 1 keV. 
In the cases with Afw = 0.7AfQ, the size of &«-dependency 
appears comparable for systems with low and high magnetic 
field strengths. 

However in systems with a more massive white dwarf, 
Afw = 1.2Af0, the emission depends more sensitively on 6**, 
for any field strength _B*. In the cases with the shock at a 
relatively low radius (e.g. Afw = 1.2Af0, m = 5 g cm~^ s~^) 
the spectra are almost indistinguishable for 9, = 0°,10°; 



the = 30° cases differ from 6*, = 0° by < 1% in the 
~ 1 keV lines. At S« = 60° the enhancement of line emission 
is < 7% when B, = 10 MG, but half as great when B. = 
50 MG. A similar variation with B» occurs among cases 
with Afw = 1.2Af0 and m = 1 g cm~^ s~^ — the greatest 6'*- 
dependency appears in the systems with the greatest relative 
shock height, rs. 



4 DISCUSSION 

4.1 Summary for magnetic cataclysmic variables 

The incorporation of 2T effects in the post-shock accretion 
flows of mCVs generally hardens the emitted keV X-ray con- 
tiimum. A dipolar accretion funnel further hardens the spec- 
trum of a 2T flow, relative to models with simpler planar 
geometry. A 2T dipolar model predicts a harder spectrum 
than a 2T planar model, if the white dwarf is as massive as 
lAf0 or more. The differences between planar and dipolar 
2T accretion models are greatest in systems where the shock 
height is large compared to the stellar radius: i.e. greater Afw 
or lower B». The spectral predictions of IT and 2T dipo- 
lar models differ more signiflcantly when the white dwarf is 
massive, has a strong magnetic field or low specific accretion 
rate. Line emission is typically weaker around ~ 1 keV and 
~ 7 — 10 keV in a 2T fiow than in a IT fiow, except for the 
Fe Ka lines, which could become stronger. 

This is qualitatively consistent with the 2T effects found 
for basic planar accretion models with out the inclusion o f 
gravity in the hydrodynamic equations (|Saxton et al.|[2005l ) ; 
Two-temperature effects upon X-ray spectra persist despite 
the introduction of a dipolar magnetic fie ld geometry, which 
in the IT model of ICanaUe et al] (120051) soften the X-ray 
spectra relative to the results of Cropper et al.l |l99^. This 
softening due to funnel geometry is, however, smaller than 
the hardening apparent when 2T physics is incorporated. 
The 2Td model predicts the most hard photons for a given 
mass Afw, and therefore its spectral fit to an observed spec- 
trum will require a lower mass than the standard planar and 
IT models do. 



4.2 Inferred masses 

To quantify how the mass estimates are affected, we calcu- 
lated spectra for models with different Af„, with the other 
parameters fixed (6'* — 0°, m = 2 g cm"'^ s~^), for medium 
(B, = 20 MG) and strong-field (B. = 50 MG) cases. For 
a given white dwarf mass, we manually match the shape 
of the X-ray continuum by adjusting the masses in ITp , 
ITd, 2Tp and 2Td models JWul 119941 : ICanalle et al]l2005l : 
ISaxton et al]|2005l . this paper respectively). Each fit is de- 
fined by matching the ratio of spectra at 0.1 keV and 10 keV, 
with the fiattest possible ratio at intermediate energies. Fig- 
ure|9]illustrates how the 2T and dipolar funnel effects change 
spectral inferences of the white dwarf mass (Afw) from X- 
ray observations. For low-mass white dwarfs (Afw ~ O.5Af0) 
the IT and 2T spectral predictions are similar. However for 
cases with higher mass, say Afw l.OAf©, the spectra differ 
considerably (middle panels). 

For large Afw, the 2Tp models require lower masses 
than ITp models. For B, — 50 MG the ITd model requires 
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greater masses than 2Tp models; but for _B* = 20 MG the 
ITp model is more massive. The 2T dipolar models consis- 
tently need the lowest Mw (see corresponding fitted masses 
in Figure llOp . The introduction of dipolar geometry to 2- 
temperature fiows (from 2Tp to 2Td models) has hardened 
the spectra. For example, a I.OOM0 2Tp model matches a 
2Td model with Mw ~ O.QSM© (i.e. a 2% discrepancy, see 
Table [!}. The X-ray continuum of a I.OOMq ITd model is 



best matches a 2Td model with ft; O.95M0 (a 5% dis- 
crepancy, see Figure 121 lower middle panel). 

We note that the flow geometry and 2T conditions can 
affect the line spectra as well as the continuum. Here, we 
demonstrate that mass estimates obtained from the contin- 
uum fitting depends on the assumed geometry of the fiow. If 
the spectral line information were also considered, we would 
need to consider variations of the values of _B«,m,&* and 
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Figure 5. Quotients of spectra in IT and 2T dipolar models (lTd/2Td). In these cases the white-dwarf mass is 0.7 A/q and accretion 
occurs onto the pole (9, =0°). 



the hot-spot area too. For the present study, it suffices to 
show that the joint inclusion of 2T phenomena and dipo- 
lar geometry can reduce mass estimates by < 9% near the 
Chandrasekhar limit, by < 7% near a solar mass, but only 
< 1% near 0.5 M©. 

We now consider the mass estimates of iRamsavl (|2000l ) 
who determined the mass of the white dwarf in 21 mCVs 
- 8 polars and 13 intermediate polars (IPs) - using RXTE 
PCA (2-60 keV) dat a. Using the ITp model prescribed by 
I Cropper et al. I (|l999l ). iRamsavl (|2000h found a mean mass 
of O.87M0 and O.SOM© for his sample of IPs and polars re- 
spectively. Applying a correction factor based on our results 
above gives a mean mass of 0.82Mq and 0.76Mq for the IPs 
and polars respectively. 

Other affects mayb e important as well. For instance, 
ISuleimanov et al] (|2005l ) found that when they determined 
the mass of 14 IPs using RXTE PCA plus HEXTE (10- 
200keV) data the mean mass of the IPs was reduced by 
0.20Mq to O.75M0. This was explained by the fact that 
in IPs broad accretion curtains with high total absorption 
columns are thought to be present, which can affect the mod- 
elling of the continuum even at relatively hard X-ray ener- 
gies. Because polars have higher magnetic field strengths and 
hence more coUimated accretion streams, absorption may af- 
fect their X-ray continua less than for IPs. 

We note that the m ean mass of whit e dwa rfs in 
mCVs as determined by ISuleimanov et al] |200^ and 
in this paper is greater than that of isolated white 



dwarfs (0.56Mm. iBergeron et al.l ll992D . We note that 
IWickramasinghe fc Ferrariol ( 20051 ) found that white dwarfs 
with high magnetic fields have a mean mass of 0.93Mq - 
significantly higher than that in isolated white dwarfs. Our 
mass estimates remain consistent with the scenario in which 
magnetic white dwarfs are more massive than their non- 
magnetic counterparts. 



4.3 Further applications and scope 

Our 2-temperature formulation with magnetic funnel geom- 
etry is applicable to a variety of accreting stellar systems 
other than the cataclysmic variables. With substitution of 
the relevant radiation processes, the model is applicable to 
accretion onto neutron stars, and similarly it is applicable 
to field-aligned accretion onto T Tauri stars. The key as- 
sumptions are that the flow be hydrodynamic (coUisional), 
fleld-aligned, and that radiation effects can be described in 
terms of local variables. 

Our hydrodynamic formulation treats accretion flows 
that are steady and smooth: rapid temporal variability is 
omitted. Continuous flows onto some mCVs oscillate at 
~ 1 Hz, interpreted as a therma l instability of the shock , 
fsee iLanger et al.lll98ll . J,982; Chevalier fc Imamural Il982l : 



ISaxton et al. 1998). This is observed optically in only a 
mino r ity of systems (e.g. Middleditch 19821 ; [Cropper et"al] 
1986 ; llmamura fc Steiman-CameronI 19861: iLarsso 

3 Il987l , 



19891 : iRamsever et aLlll993l : lMiddleditch et aLlll997h but not 
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Figure 6. Quotients of spectra in IT and 2T models {lTd/2Td), as in Figure[5]but with A/„ = 1.0 Mq. 



yet in X-rays (e.g. IWolff eFaLlllOggl : [Christianll2000l ). Thus 
fast oscillations may not affect X-ray spectra at presently 
observable levels. The absence of a steady structure would 
preclude coUisional ionisation equilibrium, which affects the 
spectral line strengths, as they critically depend o n the den- 
sity structure (see e.g. iFujimoto fc Ishida f993; IWu et al] 
120011 ). Our model is more applicable to cases with slower 
temporal variability - flaring which exceeds the post-shock 
free-fall time - by calculating sequences of steady models 
with different m values but identical Mw , B* and 6* . 

The spectral effects of spatial variations are harder to 
assess. Mild inhomogeneities might be modelled by treating 
the accretion column as a distribution of streams with dif- 
fering m and cross-sections, then summing their spectra. 
By construction, we cannot model highly inhomogeneous 
flows: these would require highly resolved, time-dependent 
hydrodynamic simulations. In some asynchronous accreting 
systems, a continuous flow may fragment into dense blobs 
due to a magnetic drag (e.g. Arons fc Lealll980l : [Frank et al.l 
1 19881 : [Kin g 1993; Wynn & King 19951). However if the mag- 
netic field is strong enough then the blobs shred into a con- 
tinuous field-aligned fiow as they near to the stellar surface 
(lArons fc "Lealll980l ). If the field is weak, blobs may survive 
to strike the stellar surface individually. Such a shower is 
qualitatively different from a pressure-supported stand-off 
shock. However, emission from steady accretion shocks is 
universally observed in mCVs and is the major component, 
except at the softest energies where the emission from bal- 
listic blobs contributes in most cases. Spectral modelling of 



predominantly clumpy accretion is a challenge beyond the 
scope of this paper. 
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APPENDIX A: COORDINATE SYSTEM 

Following [Canalle et all (120051 ). we adopt a curvilinear co- 
ordinate system (u, w, (fi) in which tp is an azimuthal angle 
about the polar axis, u; is a path-length along a magnetic 
field line, and u measures along an equipotential transverse 
to the field lines, effectively choosing the latitude of the field- 
line's footprint on the stellar surface. The dipolar coordi- 
nates correspond to conventional polar coordinates accord- 
ing to 



cos 6 



and 



implying that 



2 4 , 

w r + ur - 
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Figure 9. Matching the X-ray continua of two different accretion models requires different Mw values. The top row compares 2T 
planar and dipolar models. The bottom row compares IT dipolar and 2T dipolar models. Masses were adjusted to achieve approximately 
proportional continua (flat quotient profile). We have fixed S* = 50MG, m = 2 g cm~-^ s~^ and S» =0°. For low (left panels) the 
inferred masses are similar; but for more massive cases the model masses differ appreciably. 
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Table 2. Sets of M^, values that give comparable spectral con- 
tinua, calculated for the ITp, ITd, 2Tp and 2Td models. The 
latter imply similar or lower masses. For these mass fits, we fix 
m = 2 g cm-2 s-i and = 0°. 
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We scale the coordinate axes so that the stellar surface oc- 
curs at r = 1. The transformation to Cartesian coordinates 
obeys 
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V ur '^ cos Lp 
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\j ur^ sin if 
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The orthogonal unit vectors, {u,w,ip), vary spatially. Lo- 
cally they transform into their Cartesian counterparts {i, j, 
k) according to 



(A5) 
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where (see e.g. lArfken fc Wegerl I2001I) the transformation 
matrix is 

1 92; I dy 1 dz 
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J_ da; J_ 1 dz 

hi '5u hi ml hi Ihi, 

1 dx 1 dy I dz 

Ti2 dw Ti2 dw Ti2 dw 
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The coefficients hi, /12 and h^ comprising the metric of the 
curvilinear coordinates system, are defined by 



h,i 



[1,2,3] 



dx 



dlujWjip] J \d[u,w,(fi] J KdlUjWjip] 

with explicit functional forms given in lCanalle et al.l ()2005h . 
For a small fluid element following a magnetic field line, its 
transverse area varies as hihs in response to the converge 



of the field. The rate of convergence of the field lines, and 
thus the lateral compression of field-channelled inflow, can 
be expressed in terms of a function 

d 



H 



dw 



ln(/ii/i3 



(A8) 



APPENDIX B: NUMERICAL SOLUTION 

We find that it is convenient to express the post-shock struc- 
ture in terms of a specific entropy variable. 



s = Pp-'' = 



hih-i 



7-1 



which has a gradient equation, 
(7 - l)As 



ds 
dw 



-(7-1 



, 7-1 



An 



7-5/2 



(Bl) 



(B2) 



Since the expression (|B2|| never changes sign, the spatial 
variation of s is guaranteed to be monotonic throughout 
the post-shock flow. The maximum value of s occurs at the 
shock; the specific entropy falls to zero at the stellar sur- 
face. Thus s is a convenient replacement for w in the role 
of independent, integration variable. Differential equations 
in terms of ds are obtained by multiplying dw/ds into (|18p . 
(|19[) . (|20p and (|2ip and then simplifying algebraic terms to 
reduce numerical round-off errors. 

The strong shock boundary conditions imply that 

^ (B3) 



and 



(B4) 



(|Canalle et al.ll2005l l. where we use a unit convention based 
on C = 1, the ste llar radius = 1 and superficial escape 
velocity = yJ'2(TM^^jR^ — 1. The shock values of s 
and Po follow from l|Bip and (|14p . with <js treated as a free 
parameter. Assuming trial estimates of the shock location, 
we integrate inwards to s ^ and check for consistency 
between this inner radius and the stellar surface, r = 1. A 
numerical root-finder iterates to a value of i?s that matches 
both boundary conditions. 

This paper has been typeset from a TpTjX/ ET^jX file prepared 
by the author. 



